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The OPERA collaboration has recently reported that neutrinos travel faster than light. We review 
the theoretical situation of constraints on violations of Lorentz invariance, focusing in particular on 
the compatibility between the OPERA results with both previous constraints and ones derived after 
OPERA. We generalize to higher order operators the recent constraint provided by the absence of 
neutrino energy loss, via electron-positron pair production at OPERA energies, and show that 
simple modified dispersion relations in vacuo within an effective field theory framework cannot 
explain OPERA measurements. 

PACS numbers: 



I. INTRODUCTION 



Lorentz symmetry is one of the grounding pillars of special relativity and hence our present understanding of the 
physical world. As such, its validity has been checked by numerous experiments over the last century. Despite no 
experimental indications of any deviations, some approaches to quantum gravity suggest that due to the fundamental 
structure of spacetime Lorentz symmetry might not be an exact symmetry of nature but merely an approximate or 
emergent symmetry at low energies that is strongly broken above some energy scale. Searches for any violation of 
Lorentz symmetry (LIV) have therefore received more and more attention over the last few years as possible signals 
of physics beyond the standard model. 

Independent from any theoretical motivation, experimental constraints on LIV are very tight and so any experi- 
mental indication of a deviation from Lorentz invariance is of high importance. Recently, the OPERA collaboration 
reported pQ a measurement of the time of arrival of non-oscillated (with a contamination of v^, v e and v e that 
is estimated to be less than 5%) over the path length ~ 730 km from CERN to Gran Sasso with average energy of 
~ 17 GeV. The arrived earlier than light would, by an amount 

At = 60.7 ± 6.9 (stat.) ± 7.4 (sys.) . (1) 

No hints of energy or seasonal/sidereal dependence were found over the energy range of the experiment. This mea- 
surement can be converted into a measure of the propagation speed of the i^'s, which yields a propagation velocity 
for non-oscillated [T] of 

Ac 

— = (2.48 ± 0.28 (stat.) ± 0.30 (sys.)) x 10~ 5 , (2) 
c 

where c is the low energy speed of light in vacuum, and Ac = v v — c. 1 Previously, the MINOS Collaboration reported 
[5] Ac/c < (5.1 ± 2.9) x 10~ 5 for 3 GeV compatible with within 2a. 

Without additional new physics, the OPERA result is a statement about the behavior of (almost) free propagating 
neutrinos. The relevant modification to the neutrino Lagrangian sensitive to these results is a modification to the free 
field kinetic term, or equivalcntly the particle's dispersion relation. Below we parametrize, using some simplifying 
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1 A few months after the publication of this paper, the Opera collaboration found two possible sources of error in their experiment. The 
two new effects shift the time of flight result in opposite directions and at this point there is no new data release by Opera. However, 
the Icarus experiment reported recently |] on a re-analysis of the most recent data taken by Opera. Icarus finds a neutrino propagation 
velocity compatible with speed of light within the experimental uncertainties. Since our work does not attempt to explain the Opera 
claim but merely elucidates the difficulties the result has with conservative modifications of the Standard Model we shall keep the 
discussion basically unchanged in what follows. 
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assumptions and previous experimental results, a commonly discussed class of LIV effects on free field propagation 
within an effective field theory (EFT) framework and explain where each of the possible modifications runs into ex- 
perimental constraints. The advantage of this procedure, with respect to other purely kinematical/phenomcnological 
approaches, is that the structure of the SM is preserved, and hence observational consequences are completely calcu- 
lable. Furthermore, the simplified EFT framework here adopted allows a transparent physical interpretation of the 
constraints casted on the salient parameters characterizing the departure from Lorentz invariance. 

Since no sidereal variations have been measured for OPERA (or any other experiment) there is no experimental 
reason we cannot assume that the dispersion relations for neutrinos are rotationally symmetric. We do this for 
simplicity, although a more complicated LIV pattern could be chosen. Such an assumption picks out a preferred 
frame, and the natural preferred frame is that of the cosmic microwave background. Since the boost factor of an 
Earth centered frame is only ~ 10~ 3 [4 with respect to the CMB frame, any rotation breaking effects generated by 
the Earth's relative motion with respect to the CMB will be naturally suppressed by a factor of ~ ICP 3 relative to 
rotationally invariant effects. 

In short, the modified dispersion relation that we consider for a particle of species i in the following is 

E 2 =P 2 +^ + ^J^- 2 , (3) 

with E and p being the energy and the momentum of the particle respectively, m the particle mass, £j = ± and M{ the 
energy scale at which LIV appears. It is, of course, possible that the eigenstates of the LIV terms do not correspond to 
the energy eigenstates for neutrinos. We will not analyze the OPERA results of such a case in detail here since, as we 
shall see, existing neutrino constraints are so tight that all neutrino states must have essentially identical dispersion 
relations thereby rendering flavor dependent dispersion an unnecessary complication. The dispersion relation Eq. ^ 
has also been often considered in the literature as a test model for how quantum gravitational effects might influence 
infrared physics [SHE] . For an exhaustive classification of neutrino operators that can be added to the standard model 
see [TO]. 

We leave the index n as a free phenomenological parameter and consider the cases n = 2,3,4. The case n = 1 
would produce huge effects at low energy and is strongly constrained. We shall show that due to the strong energy 
dependence n = 4 is in conflict with the OPERA results and so any n > 4 is also ruled out. LIV enters only through 
CPT odd operators at n = 3, while CPT invariance is preserved if n = 2 or n = 4 (while CPT breaking implies 
Lorentz breaking in any local EFT - see [IT] - the reverse is not true and Lorenz violation does not imply CPT 
violation). Therefore coefficients for particles and anti-particles are related by £j = (— [15] . 

One may wonder if more general dispersion relations, in which for example various terms with different n are 
considered to contribute together, can be studied. However, this would require a strong fine tuning. Indeed, if the 
violation of special relativity comes in from new physics at some energy M, then the natural contribution of each 
Lorentz violating correction term changes the speed of the neutrino by a factor c n p n ~ 2 /M n ~ 2 , where n is the order of 
the dispersion correction and c„ is a dimensionless coefficient. If M is not near OPERA energies, then the c n 's must 
be very carefully chosen, i.e. fine tuned, to contribute equally to the change in the neutrino speed. The alternative is 
to have M, the scale of the new physics, to be at OPERA energies such that the c„ can all be of the same (perhaps 
small) size. This would, however, be problematic as physics at 10 GeV is obviously exceptionally well explored and 
one would still have to account for the smallness of the c„. 

Therefore, we assume here that there is a hierarchy of terms governed by new physics at some mass scale M and 
check each order correction in energy, which corresponds to the assumption that M is not at 10 GeV and that the 
coefficients are not fine tuned, so that naturally there is a dominant term at any given neutrino energy. Without 
any custodial symmetry, one would expect the relevant operator that generates n = 2 changes in the dispersion to 
dominate (see e.g. [T2"]). We also check n = 3 and n = 4 as several theoretical models suggest this kind of dispersion 
relations to be the dominant terms, n = 3 or n = 4 dispersion corrections would, of course, require some other 
unknown physics to prevent the n = 2 term to be dominant (see again |12j and |13l 114] for a recent discussion of these 
issues) . 

Several constraints have been cast on LIV in different particle sectors using this approach. For example, photon 
dispersion is constrained by absence of birefringence to the level M > 10 14 Mpi in the case n = 3 [TBI EL an d a 
constraint of order M > 10 4 Mpi for n = 4 is expected from ultra-high-energy 7-rays [TS]. Electron/positron dispersion 
relation is constrained by observations of synchrotron radiation from the Crab Nebula to a level M > 10 5 Mpi for 
n = 3 [19 , while again ultra-high-energy cosmic ray observations help constrain the n = 4 case to M > 10 3 Mpi [18] . 
Also in the hadronic sector constraints coming from ultra-high-energy cosmic rays are extremely strong, at the level 
of M > 10 10 A/pi (n = 3, [2D]) or M > 10 3 (n = 4 [3TJ) for protons. In the case n = 2 there no energy scale to compare 
to and constraints are directly on the dimensionless parameter Ac/c. They are numerically very strong, for example 
Ac/c 7e < 10" 13 [221 [23]. 
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At these orders in LIV dispersion, neutrinos also are tightly constrained. In this note we review previous constraints 
on LIV in the neutrino sector, discuss new ones, and detail where the recent OPERA measurement runs into direct 
confrontation with other constraints. 



II. RELEVANT NEUTRINO PROCESSES 



The relevant constraints on LIV in neutrino physics for OPERA come from flavor oscillations, time of flight ob- 
servations, and ultra-high energy neutrino observatories. In this section, we give some background on each process 
before moving on to constraints. 



A. Flavor oscillations 



Neutrino flavor oscillations yield constraints on LIV differences within the neutrino sector. Neutrino oscillations 
depend on the differences in E — p between different neutrino eigenstates. In standard neutrino oscillations, this 
difference is governed by the squared mass differences between the energy eigenstates. With LIV (and our assumption 
that the LIV eigenstates are the energy eigenstates) oscillations are governed by the differences in the effective mass 
squared, denoted N? = mf + ^p n /M" -2 . Therefore, neutrino oscillations do not probe any absolute LIV in the 
neutrino sector, but rather the differences in any LIV dispersion relations between different neutrino states. 

Let us consider a neutrino produced via a particle reaction in a definite flavor eigenstate / with momentum p. 
The amplitude for this neutrino to be in a particular energy eigenstate i is represented by the matrix Un, where 

UjjUg = Sjj. The amplitude for the neutrino to be observed in another flavor eigenstate J at some distance L 
from the source, after some time T is then 

Mj = £ Uy-^-^Un « U^e-^^Un . (4) 

i I 

The transition probability can then be written as 

^ „ / 5N^L\ , ( SNlL\ 



i,j>i 



with SNfj = Nf — Nj and fjjij and Gijij are functions of the mixing matrixes. In the standard formalism used by 
experimentalists, 



a \ LIV 
Ac 



c 



5Nl =Aml+p 2 [^) . (6) 



where now 



A \ LIV -, / \ n-2 , / \ n-2 

Ac \ n — 1 I p \ n ~ \ I p 



A nice summary of neutrino oscillation observations, with particular attention to LIV, can be found in [24] . 



1. Solar neutrinos 



Neutrinos produced by the Sun at ~MeV energies yielded the first hint of neutrino oscillations. In the LI framework, 
their flux can be understood after accounting for oscillations with Ar/jg ~ 7.58 x 10~ 5 eV 2 [551 In principle solar 
neutrinos can also be used to constrain LIV effects. However, in this case the Mikheev-Smirnov-Wolfenstein (MSN) 
effect [571 12H] must also be taken into account. Since other constraints prove more difficult to reconcile with the 
OPERA result, we shall ignore solar neutrinos, although in principle the analysis could be done. 
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2. Atmospheric neutrinos 



Muon neutrinos and antineutrinos are produced in interactions of cosmic rays with the Earth atmosphere. Experi- 
ments detect preferentially and D^. In this case, the survival probability is 

2, n x 2 / AmlL f Ac\ LIV El\ 

(2 « 3m (t + (vL^)' < 8 » 

where 0^ are the mixing angles of the neutrino mixing matrix. Best fit values (without LIV) are: sin 2 (#23) = 0.42 
and Am atm ~ 2.35 x 10" 3 eV 2 gSJUS]. 



3. Reactor antineutrinos 



Electronic antineutrinos produced by nuclear reactors with ~MeV energy also provide relevant oscillation measure- 
ments. 

Long baseline The survival probability for long baseline experiments (e.g. KamLAND [53], with a baseline of about 
180 km) is 



/ Am%,L 




LIV 




[ AE + 




long 





~ 1 - sin 2 (2fl 12 ) sin 2 + \— — 1 . (9) 

Best fit values obtained in the standard LI framework are: sin 2 (#i2) ~ 0.3 and Anig = 7.58 x 10 -5 eV 2 [2"S"Il26| . 

Short baseline Evidence of electron antineutrino disappearance was sought on much shorter baseline (L < 1 km). 
In this case the survival probability is ruled by a different set of parameters |26j 

■(-)*'(^-(t)1¥)- 

Experiments reported so far null results so far (i.e. no disappearance), indicating #13 <C 1 (in facts, hints of 
sin 2 (0i 3 ) ~ 0.02 were discussed in [25]). This also means that this observable is not sensitive (yet) to modified 
kinematics, because the small amplitude of the oscillation is controlled by #13 -C 1. 



4- Accelerator neutrinos 



At energy > 1 GeV a few experiments with short baselines L ~ 1 km have provided evidence for various oscillations, 
including —¥ v e , —¥ v T , v e — > v T and their conjugates. The T2K Collaboration, with longer baseline, reported 
recently evidence for oscillation — > v e [30 at ~ 600 MeV. This process is controlled, as is the short baseline reactor 
case, by the small angle 6*13 [26], hence it cannot be used to cast constraints on (Ac/c) LIV . Also MiniBooNE reported 
the detection of v e . Interestingly MiniBooNE finds a 3er excess at 300-500 MeV of v„ — > v l: but only a 1.3a excess in 
the conjugate channel [31 , hinting perhaps at CPT violation. MiniBooNE also searched for sidereal dependence of 
the v e signal, placing strong constraints on some combination of SME parameters |32j . 



5. Constraints 



The best constraint to date comes from survival of atmospheric muon neutrinos observed by the former IceCube 
detector AMANDA-II in the energy range 100 GeV to 10 TeV [33], which searched for a generic LIV in the neutrino 
sector [31] and achieved (Ac/c)y < 2.8 x 10~ 27 at 90% confidence level assuming maximal mixing for some of the 
combinations i,j. Given that IceCube does not distinguish neutrinos from antineutrinos, the same constraint applies 
to the corresponding antiparticles. The IceCube detector is expected to improve this constraint to (Ac/c)y < 9x 10 -28 
in the next few years [35] . The lack of sidereal variations in the atmospheric neutrino flux also yields comparable 
constraints on some combinations of SME parameters [36] , 
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B. Time of flight 



We caution the reader before starting this section: From here on out we attribute a definite velocity to the neutrino 
flavor eigenstates, although they are not energy eigenstates (or even perhaps LIV eigenstates). However, given the 
above constraints on the LIV differences between neutrino energy eigenstates, and given that we are considering 
ultra-relativistic neutrinos, for which the mass term has negligible effect over the CERN to OPERA distance, we can 
safely refer to the velocity of a flavor eigenstate. 

The dispersion relation Q implies, assuming Hamiltonian dynamics, that the propagation speed of a particle is 

TO 2 ^ ( ft — \ \ ( P \n-2 

w (p)~l-_+£ _- (-£) . (11) 



2p 2 \ 2 J \M 

In turn the time delay (or advance) upon arrival over a path length L, with respect to a light ray traveling at c = 1, 
is 

. ™, s v(p) — c L . 
AT(p) = . (12) 



We can define 

Ac v(p) — c m 



C (n) c 2p 2 2 

Given that the standard convention is to include all the factors (n — l)/2 into M, we do the same in the following 



formulae, that can then be derived from Eq. (13) by replacing M n ~ 2 — > M n ~ 2 /(n — l)/2. For n = 2, Ac/c is a 
parameter entering directly the modified dispersion relation. Given that Ac/c is the quantity "directly" accessible to 
the experiments at a given energy, we relate the mass scale M to Ac/c as 

/Ac™- m 2 V V(n_2) 

M=px hr (n) + v ' (14) 



where the sign factor £ has been absorbed in the definition of Ac/c to be positive. It is clear from Eq. ( 12 1 and Eq. ( 14 ) 
that the constraint placed on M by the measurement of a time delay depends on both the energy and propagation 
distance. 



1. Constraints 



Unfortunately we have to date only a single astrophysical event for which TOF constraints can be effectively cast on 
(Ac/c) TOF . The explosion of SN1987a was a peculiar event which allowed to detect the almost simultaneous (within 
a few hours) arrival of electronic antineutrinos and photons. Although only a handful of electronic antineutrinos at 
MeV energies was detected by the experiments Kamiokall, 1MB and Baksan, it was enough to establish a constraint 
(Ac/c) TOF < 10 -8 [37] or (Ac/c) TOF < 2x 10 -9 [35] by looking at the difference in arrival time between antineutrinos 
and optical photons over a baseline distance of 1.5 X 10 5 ly. Further analyses of the time structure of the neutrino 
signal, in particular using the fact that the least energetic neutrino in the signal (at 7.5 MeV) was detected within 10 
s from the most energetic one (at 30 MeV), strengthened this constraint down to ~ 10 -10 [3"5l 140) . 

The scarcity of the detected neutrino did not allow the reconstruction of the full energy spectrum and of its time 
evolution. In SN models two main effects are present and lead to the final time-energy structure of the spectrum. 
First, the SN is globally cooling, hence the average neutrino energy decreases with time; on the other hand, neutrino 
diffusion in the SN medium depends strongly on the energy of the neutrino, which then determines its escape time. 
Given these uncertainties, we find constraints purely based on the difference in the arrival time with respect to photons 
more conservative and robust. 

Hence, adopting Ac/c < 10~ 8 , the SN constraint implies M > 10 8 x 10 MeV ~ 10~ 13 M Pl for n = 3 or M v > 
0.1 TeV ~ 10~ 17 Mpi for n = 4. 



C. Pure LIV effects 



LIV introduces, within an EFT description, new processes that a neutrino can undergo. Since neutrinos couple 
to gravity and photons (via their magnetic moment and charge radius couplings), superluminal neutrinos will emit 
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graviton and photon Cerenkov radiation in vacuum. In addition, high energy superluminal neutrinos will emit 
neutrino/antineutrino pairs via a neutral current interaction if n > 2. The detailed rate computation for these 
reactions can be found in [30] |3T] , we merely summarize below. 



Cerenkov radiation in vacuum (photon emission): v — > vy. This possibility has been already investigated for 
renormalizable operators [32] . Although the rate of this reaction has always been considered too small to produce 
significant effects even on cosmogenic neutrinos at 10 20 eV, the effects implied by OPERA are at much lower 
scales than Mp\. Therefore, the rate can be strongly enhanced. The energy loss rate was computed in [20] and 
found to be 



M\ 2 f E \-( 2 " +1 ) 



Mpi J V 1 PeV 



10 2tm -* b s . (15) 



Gravitational Cerenkov radiation can in principle be considered. However, it is sublcading with respect to 
photon emission and can possibly be an effective energy loss process only for GZK neutrinos with E > 10 19 eV 
[46) . Furthermore, for n > 2 superluminal neutrinos, neutrino splitting dominates both Cerenkov reactions. 

Neutrino splitting (v — > vvv) was studied in the context of ultra- high-energy cosmogenic neutrinos [41] . While 
the calculation was done for n — 4, the methodology is readily adapted to any n > 2. In a LI scenario the energy 
threshold for ^-splitting would be infinite. However, with LIV there is instead a finite energy above which this 
reaction can happen. The threshold equation for this reaction is analytically solvable when all the neutrinos 
involved are in the same mass eigenstate, for which the rate is maximal. In this case the energy threshold 
goes as E t h = (m 2 ,M n ~ 2 ) 1 / n . If the outgoing neutrino/anti-neutrino pair are in a different eigenstate, then the 
threshold will change slightly but can still be solved for numerically. 

Using the relation £ p = (— l) n £„, the typical energy loss time scale for a high energy neutrino with energy well 
above the threshold energy scales as 2 

64tt 3 /A^ 3( "- 2) 



T.-splitting ^ 3G 2^ 5 { E ) , (16) 

where Gp is the Fermi constant. For example, for the lowest n for which this reaction is allowed, n = 3, this 
corresponds to 

„io-w 3 m" s . d7) 

\M P J \lOGeVj v ' 

We caution the reader that there are 0(100) phase space factors which vary for each n, so that the lifetime is 
only approximate. However, since the scaling with M is so strong, these phase factors are largely irrelevant as 
they change M by only an 0(1) factor. 

Neutrino pair production (v — > ve + e~) has been recently proposed in [33]. The calculation was focused on the 
case n = 2 as 5 was considered constant when computing the terminal energy for Opera, although much of the 
calculation and result extends to higher n with only slight modifications. In particular, the rate for neutrino 
pair production of unmodified electrons will scale the same way with energy as the rate for neutrino splitting 
far above threshold. The threshold energy will in general depend on n, and can be recovered from a scaling 
S — s- £ V (E/M) n ~ 2 only up to 0(1) numerical factors. 

The threshold equation reads (for an electron-positron pair of opposite helicity) 



™ 2 ^ +w KH)' (i8) 

where x,y and t are the fraction of initial momentum p carried respectively by the outgoing v, by e~, and by 
e + ,i=l — x — y and < x,y, t < 1. 



M r 



2 As pointed out by Ward 1421 . a dimensionless factor of (E/A4) n 2 was missing in the original rate computed in |41| . 
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The computation of the threshold is straightforward but the general solution is quite cumbersome. Henceforth, 
we shall provide here only the form of the threshold energy for n = 3 and furthermore assume m„ ss 0. We 
also assume no violation in the electron sector £ e — 0. This is justified because using only the constraints that 
can be cast on LIV in the electron-positron sector at OPERA energies, any LIV in the electron sector must be 
much smaller that the one reported by OPERA, at least for n — 2,3. 

With these assumptions, the electron/positron pair takes most of the total momentum, so that x ~ at 
threshold. This implies that the form of the threshold does not depend on n, apart from the scaling of the term 
p n /M n ~ 2 . As a result the following general scaling holds 

E th,M = ~T^~ ' ( 19 ) 

with the replacement #(„) = £ v {E th /M) n ~ 2 . 

In addition the rate of this reaction as computed in |44j is general for any n up to numerical factors once 
one performs the same substitution as before with generic energy 8 — > £, V (E / M) n ~ 2 . The generic energy loss 
time-scale then reads 

TV— pair — G"p E f J , (20) 



where we dropped the purely numerical factors. As we see, the rate matches the neutrino splitting rate Eq. ( 16 ) 
up to numerical factors when n > 2. The main difference between the two reactions is that pair production is 
allowed when n = 2 while neutrino splitting is expressly forbidden for flavor blind Lorentz breaking. 

Pair production alone yields strong constraints on OPERA and so we neglect neutrino splitting as it can only 
mildly strengthen any constraint. By integrating the energy loss rate from pair production over the distance 
CERN-Gran Sasso L we obtain 



£-3„ + l _ E - 3 n + l = (3n _ 1)eE -^-2) k j3_ L ^ E -Sn + 1 {n) 

where E is the energy on a neutrino starting with energy Eq after propagation over the distance L and E Te { is 
the energy at which we normalize the parameter £„. For Opera, £„ ~ 5 x 1CP 5 for E Te f ~ 10 to 30 GeV. We do 
this in order to have a direct connection with the observed values. The factor k = 25/448 was computed in |44j . 
The "termination" energy Et corresponds to the energy that a neutrino would approach after propagation over 
a distance L if it started with Eq ^> Et- We remark here that the termination energy Et is a mildly varying 
function of n and of the energy scale E xe f. 

Whereas the oscillation and time of flight supernova constraints are reasonably self-contained, the decay rate 
constraints are much more tightly tied to the exact characteristics of OPERA and so we devote a larger section to 
the combined constraints. 



III. CAN WE RECONCILE OPERA WITH CURRENT CONSTRAINTS? 



A. Energy dependent LIV 

It is clear that an obvious problem is to reconcile the OPERA result with the SN constraints, as they are both 
TOF limits, so they effectively constrain the energy scale at which possible LIV might be found. Ac/c < 10 -8 at 
~10 MeV is incompatible with Ac/c ~ 10 -5 at ~10 GeV if LIV is energy independent, as this would be a CPT 
invariant scenario and LIV differences between various states are so tightly constrained from oscillation data. This 
alone completely rules out n = 2 modified dispersion relation if it has to be valid at all energies, and has additionally 
been clearly ruled out at OPERA energies alone by [44j . 

On the other hand, we are also forced to impose n < 4. The reason is the following: the OPERA collaboration 
also measured the neutrino propagation time by dividing the sample of "internal" events, for which it is possible to 
measure the neutrino energy, in two subsamples with different average energies (Ei) — 13.9 GeV, (E2) — 42.9 GeV. No 
statistically significant differences were found in the time delays of the two samples: Ati = 53.1±18.8 (stat.)±7.4 (sys.), 




E (GeV) 

FIG. 1: Interpretation of time delays with energy dependent LIV. Blue: n — 4 LIV, with M = 3.5 TeV. Red: n = 3 LIV. Solid: 
lines drawn to reproduce the time advance obtained with the full sample of average energy 17 GeV (and with M = 700 TeV 
for the n = 3 case. Dashed: attempt of a fit with n — 3 that might reproduce also the measurement at the 2 different energies, 
with M = 10 3 TeV. In the n — 4 case the fit is basically impossible. Errors on the x-axis are inferred from the characteristics 
of the sample, not given by the experimental Collaboration. 



At 2 = 67.1 ± 18.2 (stat.) ± 7.4 (sys.). If we now compute the time delay on a baseline of ~ 730 km assuming n = 4 we 
find time delays incompatible with the observations, as shown in Fig. Ill 

Therefore, we are left only with n = 3, which, according to Fig.^1 can possibly explain all the time of flight 
constraints assuming M ~ 10~ 13 Mpi. The agreement with OPERA data is however achieved only at the 2a level. 
A better agreement might be achieved if we could convolve the energy dependent TOF prediction with the actual 
OPERA neutrino spectrum. Notice that, as we mentioned in the Introduction, the case n — 3 includes CPT odd 
LIV, which implies = — £ p . Therefore superluminal neutrinos imply subluminal antineutrinos. If moreover we 
assume flavor independence, we do not alter flavor oscillations. In the n = 3 case, the absence of significant energy 
dependence is interpreted as due to poor statistics. It is an easy exercise to verify that in the n = 3 case the SN 
constraint Ac/c < 10 -8 is respected. 

However, superluminal neutrinos can possibly emit e + e~ pairs, or they can split according to v — > vvv (since anti- 
neutrinos are necessarily subluminal), thereby losing energy. According to Eq. (21 1, in the case n = 3 the argument 
by [H] applies albeit with Eth — 1-5 GeV and terminal energy Et — 15 GeV, which again rules out at OPERA 
energies the most mildly energy dependent dispersion relation n = 3. 

Furthermore, according to Eq. (16), atmospheric muon neutrinos at 10 TeV would h ave a typical mean free path 



of about 3 x 10 cm (the mean free path implied by Cerenkov is 3 x 10 cm, see Eq. 15 1 for M ~ 10 13 Mpj, too 



short to be compatible with the observation of upward-going atmospheric neutrinos up to 400 TeV by IceCube. 3 
Therefore, LIV in EFT with any energy dependence does not seem a viable scenario. 



B. Flavor dependent LIV 



Given that the flavor independent LIV scenario has problems, one might be tempted to increase the degrees of 
freedom of the problem and reconsider a form of LIV that depends on flavor. In this way, the SN constraint does not 
need be necessarily respected by OPERA neutrinos, because it is cast on v e , while OPERA observes v^. However, 
the constraints on n = 2 |44j and on n — 3 just seen in the previous section still apply also in the flavor dependent 
case. Hence flavor dependence cannot be considered a solution. 

Nonetheless, we want to briefly show here that even forgetting the above mentioned constraints, it is still possible 
to rule out both n — 2 and n — 3 possibilities simply on the basis of neutrino flavor oscillations. 

First of all, given the atmospheric neutrino constraints reported in section (II A 5), we have to demand that the 



3 While submitting this paper we found that a similar remark was also made in [45] . 
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FIG. 2: Survival probability P^ e ,^ e for different values of (Ac/c) LIV . The red solid curve is the standard LI case with solar 
neutrino parameters. Green dotted and black dashed curves are for {Ac/c) LIV = 10" 16 and 5 x 10" 16 respectively. 



LIV parameters be closer than 10~ 27 for n — 2, or 10~ 31 for n = 3 (see Eq. ^ if M ~ 10~ 13 M P \ for ~ v T and 
— v r oscillations. 

We now consider reactor antineutrinos. For our purpose we can use Eq. (|9| to understand what happens in the 
LIV case. In Fig. 2 we show how even a tiny modification Ac/c ~ 10~ 16 would lead to large unacceptable departures 
from the LI case. Indeed, KamLAND observed [29] the first 2 humps of the survival probability at values L/E ~ 30 
to 50 km/MeV, in clear disagreement with the 4 humps present for Ac/c ~ 10~ 16 . Given the TOF constraints, in 
the computation of LIV effects on oscillations we have to require (Ac/c)^ ~ 10~ 5 for n = 2 and (Ac/c)^ ~ 10~ 8 
for the CPT odd n — 3 scenario at MeV energies. Hence, these values are far too large to be compatible with the 
observed oscillations. Therefore a flavor dependent LIV case is strongly disfavored. 

IV. CONCLUSIONS 

In this note we discussed the OPERA results. Our conclusions, in agreement and complementing other similar 
analyses [HI U7J HB]> are that these results are not compatible with other observations and derivable constraints 
within an EFT framework (even models predicting modified dispersion relations on finite energy ranges, including 
OPERA energies [48]). 4 

Of course, an experiment is not proved wrong by the theoreticians' inability to come up with viable explanations, 
and the last word on the reliability of the OPERA measurement lays in the hands of forthcoming experiments aimed 
at reproducing the results. However, we think that taken at face value, this analysis (and the above mentioned similar 
ones) should be taken as a strong case for being cautious in developing freely new exotic models aimed at fitting 
OPERA data, or at least as guidance about what these models should actually provide. 

It is difficult to envisage which models can escape the wealth of constraints at these energies, in particular the 
one provided by the absence of pair production by OPERA neutrinos. One example could be models like Deformed 
Special Relativity (DSR) which do not lead to new threshold reactions [53 . 5 

Alternatively, one can envisage the possibility that in vacuo neutrino propagation is not affected itself by LIV, while 
some LIV interaction can be the source of OPERA anomaly. In the past there have been models, most noticeably 
the so called space-time foam ones, like [55], predicting a modified propagation only in some SM sector (photons and 
neutrinos in [55] ) and only on average, not instantaneously, due to many interactions with some hidden sector of the 
theory (in the case of [55] . string theory D-branes). Averaging over many of these interaction events leads to features 
not predictable via standard EFT reasoning. However the model described in |55j predicts subluminal propagation of 




4 See also |49l I50| for alternative albeit less stringent arguments. For tentative explanations see e.g. |51l I52j . 

5 Unfortunately, the various DSR candidates face in general major problems regarding their physical interpretation and a working model 
is not yet available (see however [54_. for recent attempts in new directions). 
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photons and neutrinos 6 , and hence it is not amenable to explain Opera results, although it is possible that different 
implementations of the idea that lead to superluminal neutrinos may yield a workable proposal. 

Models of this kind could possibly evade most if not all the constraints based on particle interactions (see however 
|56j for a constraint on a possible incarnation of the framework discussed in [55] ) and still lead to modified propagation 
on sufficiently long distances. Furthermore, the "hidden" fields can be homogeneously distributed in the Universe, or 
highly clumped, e.g. around the Earth's position. This could perhaps play an important role, like for example in the 
model of [57] . If OPERA is confirmed, there is a long road ahead for model builders trying to match all the existing 
constraints. However, precisely because of this wealth of data we have, there is a lot of experimental guidance on how 
one should proceed. 
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